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Abstract 

We compute mean field phase diagrams of two closely related interacting fermion 
models in two spatial dimensions (2D). The first is the so-called 2D t-t'-V model de- 
scribing spinless fermions on a square lattice with local hopping and density-density 
interactions. The second is the so-called 2D Luttinger model that provides an effec- 
tive description of the 2D t-t'-V model and in which parts of the fermion degrees of 
freedom are treated exactly by bosonization. In mean field theory, both models have 
a charge-density-wave (CDW) instability making them gapped at half-filling. The 2D 
t-t'-V model has a significant parameter regime away from half-filling where neither 
the CDW nor the normal state are thermodynamically stable. We show that the 2D 
Luttinger model allows to obtain more detailed information about this mixed region. 
In particular, we find in the 2D Luttinger model a partially gapped phase that, as we 
argue, can be described by an exactly solvable model. 



1 Introduction 

This is the second in a series of papers that aim to develop a method to do reliable com- 
putations in a spinless 2D lattice fermion model of Hubbard type; see [1] for a concise 
summary. In the first paper in this series [2], an effective model for the low-energy physics 
of the lattice system was derived. It was shown that parts of the fermion degrees of freedom 
in the model can be treated exactly using bosonization. In this paper, we apply mean field 
methods to the remaining degrees of freedom. We present analytical and numerical results 
showing that our method is useful for obtaining quantitative physical information about 
the lattice fermions. 



1.1 Motivation 



The difficulty to do reliable computations in 2D lattice fermion models of Hubbard-type 
has remained an outstanding challenge in theoretical physics for many years. One can hope 



1 



that a solution to this problem will be a key step towards a satisfactory theory of high- 
temperature superconductors [3]. One simple example is the so-called 2D t-t'-V model 
describing spinless fermions on a square lattice with local hopping and density-density 
interactions; see Section 12.11 for a precise definition. One of us (EL) proposed a particular 
partial continuum limit of this lattice system [1] [2]. This leads to an interacting model of 
so-called nodal fermions, which have linear band relations, coupled to so-called antinodal 
fermions with hyperbolic band relations. It was found that this is a natural 2D analogue of 
the Luttinger model, not only in that it arises as a continuum limit of the 2D t-t'-V model, 
but also since the nodal fermions can be bosonized and thus treated exactly. In particular, 
it is possible to integrate out the bosonized nodal fermions and thus obtain an effective 
model for the antinodal fermions; see Section 12.21 

In this paper we use mean field theory to address the question if and when the antinodal 
fermions in this 2D Luttinger model have a gap. These are key questions since if the 
antinodal fermions are gapped they do not contribute to the low-energy physics. We then 
obtain an effective Hamiltonian of nodal fermions that is exactly solvable. We find a 
significant parameter regime away from half- filling where this is indeed the case. As a 
motivation for our work, we also present mean fields results for the 2D t-t'-V model that 
show that there is an interesting region away from half filling where a direct application of 
mean field theory fails. We find that this regime becomes accessible to mean field theory 
by using the 2D Luttinger model. 

We recall the key parameters of these two models. The 2D t-t'-V model is characterized 
by the nearest-neighbor (nn) hopping constant t > 0, the next-nn (nnn) hopping constant 
—t/2 < t' < t/2, the nn density-density coupling strength V/2 > 0, and the filling factor 
< u < 1 (see Section [2.11 for more details). The 2D Luttinger model depends on two 
additional parameters, k and Q, which have the following significance. To derive the model, 
it is assumed that there is an underlying Fermi surface that, when the system is near half- 
filling, is a line segment in each nodal region. These are the so called nodal Fermi surface 
arcs, and the parameter < k < 1 determines the size of these arcs. Furthermore, the 
parameter Q ~ vr/2 fixes the point about which the nodal band relations are linearized; the 
details are given in Section 12.21 One important question addressed in this paper is how to 
fix the parameter Q. 

1.2 Mean field phase diagrams 

Two-dimensional lattice fermion systems with repulsive interactions are often insulators 
at half-filling, but away from half-filling there are competing tendencies that lead to rich 
Hartree-Fock (HF) phase diagrams. A well-known example is the 2D Hubbard model which, 
at half-filling, has an insulating antiferromagnetic (AF) HF ground state However, 
away from half-filling, unrestricted HF theory yields intricate solutions that include domain 
walls, vortices, polarons etc. imposed on an AF background; see e.g. [5] and references 
therein. These solutions suggest that the pure AF state is only stable at half-filling and 
that additional holes or particles tend to distribute so as to perturb the AF state as little as 
possible. Furthermore, away from half-filling, these solutions break translational invariance 
in a complicated manner and are highly degenerate. One thus expects that the low-energy 
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properties of the lattice model away from half-filling should be dominated by fluctuations 
between these degenerate solutions. Unfortunately, it is difficult to formulate a useful low- 
energy effective model for this situation. 

Unrestricted HF theory is computationally demanding and thus applicable only for 
moderate lattice sizes. However, it is possible to find "mixed" regions in the phase diagram, 
i.e. regions with intricate and highly degenerate HF solutions, already in mean field theory. 
By the latter we mean HF theory restricted to states invariant under translations by two 
sites [HI [7] . This requires little computational effort and thus allows to explore the full phase 
diagram for arbitrarily large lattice sizes. In this manner one can identify the regions in 
the phase diagram in which an exotic physical behavior can be expected. It is remarkable 
how rich the resulting mean field diagrams for the 2D Hubbard model are [7]. 

In this paper we present mean field phase diagrams for the 2D t-t'-V model, which is a 
spinless variant of the Hubbard model and therefore somewhat simpler. We find a stable 
charge- density- wave (CDW) mean field ground state at half-filling u = 0.5, a translation- 
invariant normal (N) state far away from half-filling, and extended doping regions with 
mixed phases in-between. We indicate the mixed phases in our phase diagrams by horizontal 
lines to emphasize that mean field theory fails here, i.e. it does not allow any specific 
conclusions to be drawn; see Figure [T](a). Put differently, the mixed regions are "unknown 
territories" (analogous to white spots on ancient maps) for mean field theory, but the 
knowledge of their existence is nevertheless important physical information. 

Our main result is to show that the 2D Luttinger model provides a tool to explore these 
"unknown territories" using mean field theory for the antinodal fermions. To be more 
specific, we compute mean field phase diagrams for the effective antinodal Hamiltonian 
defined in Section 12.21 We find that the antinodal fermions indeed have a CDW gap in a 
significant region of the parameter space, as conjectured in fT, We refer to this also as a 
CDW phase, but we emphasize that, in general, it corresponds to a partially gapped phase 
of the 2D Luttinger model. There are also the nodal fermions that are gapless, and these 
fermions can dope the system even if the antinodal fermions remain gapped and half-filled. 
In this way a large part of the mixed regions of the phase diagram for the 2D t-t'-V model is 
filled in; see Figure [T|^b). We also study how sensitive the occurrence of a partially gapped 
phase is to variations in the model parameters, especially k and Q. 

1.3 Related work 

The derivation of the 2D Luttinger model was inspired by important work of Mattis [S], 
Schulz [H], Luther [TU], and Furukawa et al. [TT]; see [2] for discussion and further references. 

The phase diagram of the 2D t-t'-V model, at and away from half-filling, has been 
studied using various techniques. Recent work close to ours is [12] in which the possibility 
of phase separation is investigated; see also [131 [E]. In particular. Figure 1(a) in [12j is 
the same as our Figure [T]^a). Note though, unlike [12], we do not necessarily interpret the 
horizontally lined region in our Figure [D^a) as a phase-separated state. Instead, we only 
conclude that the considered mean field theory fails to give a stable homogeneous phase 
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Figure 1: Comparison between the coupling {V) vs. filling (z/) mean field phase diagrams 
of the 2D t-t'-V model (a) and the 2D Luttinger model (b) for t = 1 and t' = 0. Shown 
are the charge-density-wave (CDW) and normal (N) phases as a function of coupling and 
filling at zero temperature. The regimes marked by horizontal lines are mixed, i.e. neither 
the CDW nor the N phase is thermodynamically stable. The CDW phase in (a) exists only 
at half-filling u = 0.5. The CDW phase in (b) corresponds to a partially gapped phase. 

there. 

One motivation for our work are experimental results on high temperature superconduc- 
tors. It is known from angle-resolved photoemission spectroscopy [15] that these materials 
can have an electronic phase in which parts of the underlying Fermi surface do not have 
gapless excitations. In particular, for hole-doped materials, one finds near half-filling that 
the degrees of freedom in the antinodal regions are gapped, while in each nodal region 
there is an ungapped Fermi surface arc; see e.g. Figure 5 in [16j or Figure 1 in [17] We 
believe that this suppression of accessible electronic states at the Fermi level seen in exper- 
iments provide strong support for our approach. However, we stress that the emphasis of 
the present work is mathematical, and a confrontation of our results with experiments is 
postponed to future work. 

^The experimental situation for electron-doped materials is different and its relation to the present work 
will not be discussed here. 
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1.4 Notation and conventions 



We use the symbol "=" to emphasize that an equation is a definition. We denote by 
Re(c) and Im(c) the real- and imaginary parts of a complex number c, and c is its complex 
conjugate. We write "ao,i = A ± 5" short for "oq = A + B and ai = A — B" etc. We 
use bold symbols for matrices, e.g. 1 is the identify matrix. By "X = 2.31(1)" we mean a 
numerical result "X = 2.31 ± 0.01". 

The fermion models considered in this paper are defined by Hamiltonians of the following 
generic type 

kl klmn 

with fermion operators c^^"* labeled by a finite number M of one-particle quantum numbers 
k. Note that N equals the number of one-particle degrees of freedom that are included in 
the model. Our normalization is such that 

{c^,4} = 5,fc. (2) 

The model parameters tki = and Vkimn = Vmnu correspond to the matrix elements of 
the kinetic energy and two-body interaction potential, respectively, and is the chemical 
potential. The model is defined on a fermion Fock space with a vacuum |0) annihilated by 
all Cfc. Expectation values with respect to a given state of the model (both zero and non-zero 
temperatures) are denoted by (■). The state can be either the exact thermal equilibrium 
state or an approximate Hartree-Fock state; it will always be clear from the context which 
is meant. The inverse temperature is denoted by /3 > 0. 

1.5 Plan of paper 

The two models we consider are defined in Section [2l Some results presented here are (mi- 
nor) generalizations of the corresponding ones in [2j , as further elaborated in Appendix |X1 
Section [3] explains the method we use to compute mean field phase diagrams, with some 
technical details deferred to appendices [B] and O Our results are given and discussed in 
Section HI Section |5] gives some closing remarks. 

2 Models 

In this section, we define the two studied models and discuss the relation between them. 
We note in passing that the notation used here is slightly different from that in [2]; see 
Appendix lA.ll for details. 



5 



2.1 2D t-t'-V model 



The 2D t-t'-V model describes spinless fermions on a square lattice with sites and lattice 
constant a = 1. The fermions hop with amplitudes t and t' between nn and nnn sites, 
and fermions on nn sites interact with a density-density interaction of strength V/2. The 
Hamiltonian is given in Fourier space by 

Hu^y = J2 W - /^]cnk)c(k) + ^ E ^k„k„k3,k4Cnki)ct(k3)c(k4)c(k2) (3) 
keBZ kjEBZ 

with the tight-binding band relation (k = {ki, ^2)) 

e(k) = -2t[cos(A;i) + cos(A;2)] - 4t' cos(A;i) cos(A;2) (4) 
and the interaction vertej|§ 

^^ki,k2,k3,k4 = M(ki - ka) 1 — k2+k3— k4,27rn (5) 

neZ2 

with 

m(p) = cos(pi) + cos(p2)- (6) 
The fermion operators c*-^''(k) are labeled by momenta k in the Brillouin zone 

BZ = {k = {ki, ^2) \ -IT <kj <7T , kj = ^{rij + i) , rij e Z , j = 1, 2}. (7) 

Other conventions used are explained in Section 11.41 (with the k there corresponding to 
our k here). Note that the number of different one-particle quantum numbers k equals 
the system volume, M = L?'. The chemical potential /i is to be determined such that the 
fermion density 

^=^,{N), iV^='5^ct(k)c(k) (8) 

keBZ 

has a fixed specified value. We refer to v as filling factor or filling, and v — 0.5 is called 
doping. The filling lies in the range < < 1, and = 0.5 corresponds to half-filling. We 
always assume —t/2 < t' < t/2, > 0, and that L/2 is an integer. To simplify notation we 
set t = 1 in all figures, i.e. energies are measured in units of t. Figure [2] shows the Brillouin 
zone of this model and a typical example of a non-interacting Fermi surface defined by 
e(k) = fi. 

Note that we use anti-periodic boundary conditions and a large- distance cutoff different 
to that used when deriving the 2D Luttinger model in [2]. This is legitimate since we are 
interested in the thermodynamic limit L — )• 00, and finite size effects are negligible for the 
system sizes {L > 100) we use in our numerical computations. 

The invariance of the 2D t-t'-V model under particle-hole transformations provides an 
important guide for us. The interested reader can find more details in Appendix lA.li 

^The sum over n in ([5]) takes into account possible umklapp processes. 
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Figure 2: Fermi surface for the tight-binding band relation in with t = 1, t' = —0.2 and 
fi = —0.51(1) (bent dashed curves). Also indicated are the eight regions labeled by (r, s) 
and discussed in the main text (large squares for r = ±, s = 0, rectangles for r = ±, s = ±, 
and small squares for r = ±, s = 2). Note that the Brillouin zone is equal to the union of 
these regions. The two antinodal points Qr.,o and four nodal points Qr,± are depicted by 
dots and rings, respectively. The parameters used in the plot are k = 0.6 and Q/tt = 0.4, 
and /i corresponds to Qo/n = 0.45(1). 

2.2 2D Luttinger model 

A detailed derivation of the 2D Luttinger model and its partial solution by bosonization 
were given in [2]. Here we first describe this model and then summarize the results from 
[2] that we actually need. 

The 2D Luttinger model involves six fermion flavors labeled by a pair of indices (r, s) 
with r = ± and s = 0, ±. These flavors correspond to different regions in the Brillouin zone 
of the 2D t-t'-V model, as shown in Figure [2] (two large squares for s = and four large 
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rectangles for s = ±). The sizes of these regions are determined by a parameter k in the 
range < /t < 1 (we used k = 0.6 in Figure [2]). There are two more fermion flavors with 
r = ± and s = 2 (two small squares in Figure |2]), but in this paper we assume that the 
parameters are such that the s = 2 fermions are far away in energy from the Fermi level 
and thus their dynamics can be ignored [2] (we discuss this point in Sectional Remark 2). 

The momenta in the regions with s = can be written as Q±,o + k with 

Q+,o = (7r,0), Q_,o = (0,7r) (9) 

and k in 

A* {k G BZ I \ki ± k2\ <KTT, h±k2 = ^(n± + I) , n± G Z}. (10) 

Close to the points Qr,o, r = ±, the band relation in (jlj) can be well approximated by the 
lowest-order non-trivial terms in a Taylor series expansion e(Q,.^o + k) ^ At' + er{\<-) with 

er{k)=rt{kf-kl)-2t'{kl + kl), (11) 

i.e. Qr,o correspond to saddle points of e(k) if —t/2 <t'< t/2. This also explains why we 
impose these bounds on t'. Similarly, the band relation for fermion degrees of the freedom 
corresponding to s = ± can be well approximated by bands linear in {ki ± ^2) close to the 
points Qr,± = (rQ, irQ) with Q another free parameter (r = + or — ). We use standard 
terminology [15] and refer to the fermions with s = as antinodal and the fermions with 
s = ± as nodal, respectively. In Figure [2] the antinodal points Q^^o are indicated by dots 
and the nodal points Qr,± by circles. 

The 2D Luttinger model is defined by a Hamiltonian of the form Hn+Ha+Hna where if„ 
and Ha include terms depending only on the nodal and antinodal fermions, respectively, 
and Hna are interaction terms with both kinds of fermions. It is obtained from the 2D 
t-t'-V model by certain approximations that amount to a particular partial continuum 
limit [2]. A key assumption is that there is an underlying Fermi surface in the nodal 
regions s = ± consisting of line segments ( "Fermi surface arcs" ) and containing the nodal 
points (rQo, ^^Qo) for some Qo ~ 7i/2 determined by /i, k and Q; see (HOll . In [2] we fixed /j, 
by the condition Qo = Q, but in the present paper we work in the grand canonical ensemble 
and thus allow /i (i.e. Qo) to be arbitrary in intermediate steps of our computations. No 
assumption is made on the Fermi surface in the antinodal regions. 

As explained in |2j. Sections 2 and 6.2, we need to restrict ourselves to parameters such 
that 

(1 - K)n/2 < Q < (1 + K)n/2 and < < 1 (12) 

27i[t + 2r cos{Q)\ 

with the same bound for Qo as for Q. Moreover, we require Qo and Q to be different from 
7r/2, since otherwise one has additional back-scattering terms in the 2D Luttinger model 
which spoil a simple treatment of the nodal fermions using bosonization; see [2]. These 
conditions define the parameter regime of interest to us. 

The nodal fermions in the 2D Luttinger model can be bosonized and integrated out 
exactly. This yields an effective model for the antinodal fermions c£''(k) = c^^^(Q±_o + k) 
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that, in the local-time approximation [2], is given by a Hamiltonian of the form 

2^ 
L2 



Ha = ^ ^[er(k) -/ia]4(k)a,(k) + ^ ^ 4i-k2+k3-k4,ocUki)cl(k3)c_(k4)c+(k2). 

keA; r=± 



(13) 

This Hamiltonian is also of the generic type discussed in Section 11.41 (with the one-particle 
quantum numbers k used there to be identified with (r, k) here) . Note that there are 
(kL)2/2 different momenta k, and that the number of one-particle degrees of freedom in 
the model is TV = (kL)^. The filling factor of the antinodal fermions is therefore 



dof 
I' a = 



^,{Na), iV/=^ J] J]ct(k)c,(k), (14) 



keA; r=± 

while the total filling factor of the 2D Luttinger model is (including the nodal fermions etc.) 



l + il-K)C-^-l) + K'iu^-ly, (15) 



see [2]. It is 1/ in f[T^ that is to be identified with the filling factor of the 2D t-t'-V model. 

A main result in |2j are explicit formulas for the parameters /Iq and Qa in terms of the 
other model parameters: 

/ia = /i - 2VZ/ - At' + gal^aK.^ (16) 

and 

ga = 2V-g,s (17) 

with 



sin(Q)7r[t + 2t' cos(Q) + ^(1 - k) sin(g)] 
The parameter Qq can be computed from the following identity 



H = V2vf{Qo -Q) + 2Vu - 4tcos(g) - 4t'cos^(g) - 2V^Ccos(g) (19) 

with 

vf = 2y2sin(g)[t + 2t'cos(g)] (20) 
and where we have introduced a convenient short-hand notation 

C (1 - n) cos(g)(2Ji - 1) + 1(1 - nf. (21) 



The constant —ges corresponds to a renormalization of the bare antinodal interaction 
2V and arises from integrating out the bosonized nodal fermions. 

The 2D t-t'-V Hamiltonian is equivalent (in a low-energy approximation) to the effective 
antinodal Hamiltonian in (1131) only if one takes into account the additive constant (see 
Appendix IA.2D 

^kin ~l~ ~l~ ^int ~l~ (22^ 
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with 



kin 



= 1(1 - /t)2(-4(t + t') + lit + 2t')(l - /t) V) + (1-k){^-1 + k) 
X {V2vf {Qo/2 - Q + tt (1 - k) /4) - 4t cos(g) - At' cos'^{Q)), 



(23) 



(24) 



and 



£int/L^ = -fii^ + Vu^ + V[k{1 - k) cos(g)]2 - VC^. 



(25) 



The constant Sn is derived in [2J, but in the present paper we only need that it is independent 
of the chemical potential: 



As we show later, the parameter Q can be fixed by the self-consistency condition Q = Qq. 
Thus the effective antinodal model contains one more free parameter as compared to the 
2D t-t'-V model, namely n. 

As already mentioned, the results in |2] are restricted to the special case in which is 
explicitly fixed by the condition Qo = Q, and they are written in a slightly different form. 
The interested reader can find details about how to obtain the results given here from the 
ones in [2] in Appendix IA.2[ 

It is not essential to work with the Taylor expansion of the antinodal band relations, 
and one can equally well use the full band relations 



instead of (fTTl) . As we will discuss, the results for the band relations in ( [TT|) and (p71) agree 
quite well for smaller values of k, but for k close to one there are some quantitative differ- 
ences. Furthermore, in the derivation of the 2D Luttinger model, certain approximations 
were done on the interaction vertex of the 2D t-t'-V model; see Section 5 in [2]. It was 
argued that this had no important consequences for low-energy scattering processes. We 
note here that these approximations were only necessary for processes that include nodal 
fermions, and it would be possible to use the full interaction vertex for processes involving 
just antinodal fermions. With this, and using the full band relation fl27|) . one could derive 
a refined 2D Luttinger model for which the 2D t-t'-V model is recovered by setting k = 1. 

We finally mention that the local-time approximation in [2] was only done for simplic- 
ity, and it is possible to generalize our treatment here to take into account the full time 
dependent interaction. We hope to come back to this in the future. 



In this section, we discuss the mean field Hamiltonians used for deriving phase diagrams 
and the procedure that allows us to identify the mixed regions. More explicit details can 
be found in appendices [B] and O 




0. 



(26) 



4o)(k)=e(Q,,o + k)-e(Q,,o) 



(27) 



3 Method 
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The fermion models considered in the previous section are defined by Hamiltonians 
of the type ([1]). Conventional HF theory at zero temperature and fixed particle number 
Nq for such models amounts to considering the set of all variational states of the form 
V = cKfi)cKf2) ■ ■ ■c^{fNo)\0) with c\fj) = X]fc(/i)fc4 fj orthonormal one-particle 
statesH Using standard terminology, we refer to the states r] as Slater determinants. The 
fj are to be chosen so as to minimize the HF energy 

nuF = {v,Hv). (28) 

It is straightforward to compute an explicit formula of ^hf in terms of the one-particle 
density matrix = Ylk=iifk)iifk)j; see Appendix E 

In the present paper, we work at (mainly) small but non-zero temperatures (unless 
otherwise indicated /3 = 10^). This means that we minimize the full grand canonical 
potential (including the entropy) with respect to all HF Gibbs states; see Appendix [B] for 
a full discussion. We will however set the temperature to zero in the current section to not 
burden the presentation. 



3.1 2D t-t'-V model 

It is convenient to choose the Slater determinant rj as ground state of a reference Hamilto- 
nian 

i^HF = Ho + J2 ^(k, kOct(k)c(kO Yl ^-(k, k')ct(k)c(k') (29) 

k,k' k,k' 

with Ho the non-interacting Hamiltonian obtained from the one in ([3]) by setting = 0, and 
w(k, k') = w{k.', k) the matrix elements of the HF potential w. This allows to parametrize 
HF theory by the one-particle Hamiltonian /i^ defined in (129|) . As explained below, it is 
important to use the grand canonical ensemble, i.e. to fix the particle number by adjusting 
the chemical potential /x [7j. One thus obtains 

7(k, k') = J2 0{-e,)Mk)j;(k^ (30) 

k 

with Ck and fk the eigenvalues and corresponding orthonormal eigenvectors of and 6 the 
Heaviside function. Note that the chemical potential is included in e^. 

Unrestricted HF theory amounts to determining the HF potential w that minimizes Q^p 
under the filling constraint 

. = -^(%Z^. (31) 

O/J, 

This method is computationally demanding and thus restricted to small system sizes. 

By mean field theory we mean the restriction of HF theory to states that are invari- 
ant under translations by two sites. As explained in Appendix IC.ll this corresponds to 
considering the following restricted set of HF potentials: 

w(k, k') = (go + qi [cos(A;i) + cos{k2)]) 4,k' + A (5k,k'+Q (32) 



^To be more precise: the fj are vectors in C"^ with components {fj)k, and J2k ifi)kifj)k 
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with Q = (vr, vr) and three real variational parameters go, Qi and A. For A = this 
corresponds to a normal (N) state that is translation invariant, and for A 7^ one has 
a charge-density-wave (CDW) state for which translation invariance is broken down to 
translations by two lattice sites. The computational problem is now very easy: there 
are only three variational parameters, and f2HF can be computed analytically by Fourier 
transformation; the interested reader can find details in Appendix IC.li 




Figure 3: Hartree-Fock energies of the charge-density-wave (CDW) and normal (N) states 
vs. chemical potential /i for the 2D t-t'-V model at t = 1, t' = 0, = 4. The energy f^HF of 
the thermodynamically stable phase is given by the minimum of these curves. The constant 
slope of the CDW curve m < ^ < ^2 implies that the CDW phase can exist only at 
half- filling, and the kinks at yU = yUi^2, where the system changes from the CDW to the N 
phase, imply mixed regions where neither the CDW nor the N phase is thermodynamically 
stable. (A convenient overall constant has been added to ^hf-) 

It is important to note that, by working with the grand canonical ensemble, we not only 
can detect variational ground states given by Slater determinants r/cDW (with A 7^ 0) or 
?7n (with A = 0), but it is also possible that the state with lowest energy is mixed and of 
the form 

A|?7cDw)(^CDw| + (1 - A)|r7N)(?7N| (33) 

for some < A < 10 In general, the two Slater determinants correspond to different filling 
factors i^cDW and z/n, and the state in ( 133|1 corresponds to total filling 

z/ = Az/cDw + (1 - A)z/N- (34) 

The possibility of obtaining (!33|) can be seen by computing the HF energies for the CDW 
and N states, f2^p for X = CDW and N, as functions of the chemical potential fi; see e.g. 
Figure [3] for t = 1, t' = and V = 4. One finds two different regimes: there is an interval 

''We identify a Slater determinant rj with the state |?7)(?7|. 
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fil < fi < ^2 where f^Hp^ < ^^hf) i-^- ^^e system clearly has a CDW mean field ground 
state r^cDw; for V/t = A and t' = we find fX2 — V = V — Hi = 3.01 The second kind 
of regime is for /i < /xi and /x > /i2 where one finds f2§p^ < ^^hf? i-^- ^ ^ ground state rj-^. 
However, not all doping values can be realized in this way: according to ( ISTj) we can obtain 
the doping u from the slope of 

r]HF = min(fi™,^SF), 

but this is not continuous at /i = fij, j = 1, 2, where the CDW state changes to the N state. 
The mean field phases can be determined by the following four doping values, 



dof 



3 



, J = 1, 2, X = CDW, N (35) 



as follows. The system is in the CDW phase in the doping regime z/cdw,i < < '^cdw,2 and 
in the N phase for z/ < z/n i and u > z/n,2- However, in the regions i/n,i < i' < t'cDw.i and 
^CDW,2 < v < z^N,2 the mixed state in (|33|) . with A determined by flMl) . has lower energy 
than either of the pure Slater states. We refer to the latter as a mixed phase. We emphasize 
that mixed phases occur in large parts of the phase diagram; e.g. for t = 1, t' = and 
1/ = 4 we find i^cdw.i = i^CDW,2 = 0.5 and i/n,i = 1 — ^'n,2 = 0.30(1). 

The results described above have the following physical interpretation (we only discuss 
the regions close to /i = /i2 since the other one is similar). In the CDW phase one has the 
effective band relations (see Appendix IC.ip 



e±(k) = -4t' cos(A;i) cos(A;2) + go - ± \/(-2t + gi)2[cos(A;i) + cos(A;2)]2 + A^. (36) 

This shows that the CDW phase has a band gap 2|A|, and as long as /i is in this gap, 
changing it cannot affect doping. Thus the HF energy is a linear function of yU with slope 
'^CDW = 0.5 (half-filling) in this region. The N phase is not gapped and doping can be 
monotonically increased by increasing /x, and therefore the HF energy is a strictly concave 
function of /i. Thus, as we try to increase doping by increasing /i in the CDW phase, the 
HF energy of the N phase decreases faster than the HF energy of the CDW phase, and 
when both energies become equal at /x = /i2 the doping i/n 2 in the N phase is significantly 
larger than the doping z/cdw,2 in the CDW phase. 

A possible interpretation of the mixed state in (!33|) is a phase-separated state in which 
parts of the system are in the CDW phase and parts in the N phase [TJ [12] . We can therefore 
conclude that, for z/n 2 < i' < '^cdw,2, a phase-separated state has lower variational energy 
than any simple mean field state. However, we emphasize that the occurrence of a mixed 
phase does not necessarily imply phase separation, but it nevertheless proves that a true 
HF ground state is very different from any state that can be described by a simple mean 
field ansatz fl32|) (a true HF ground state can in principle be found by unrestricted HF 
theory). Thus mean field theory allows to determine those regions in phase space where 
non-conventional physics (not describable by mean field theory) is to be expected. 



""The accuracy of reached in our computations is actually greater than what the numerical error given 
here suggests. A similar remark applies to other numerical results given in this paper. 
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It is interesting to note that, for non-zero t', it is possible to have a CDW phase also 
away from half-filling. This can be seen by computing a plot similar to Figure [3] but with 
t' = —0.2, for example. One still finds that the CDW energy fi^F^ ^ function of is a 
straight line in most of the interval fii < fi < fi2, but as /i approaches fi2 from the left, this 
curve starts to bend so that t'cDW,2 = 0.53(1) > 0.5. There is no bending of fi^p^ close to 
fi = fii, however, and t'cDW,! = 0.5. Thus at t = 0, t' = —0.2, = 4 it is possible to dope 
the CDW state on the particle side (i.e. for u > 0.5) but not on the hole side {i> < 0.5). 
This shows that the parameter t' affects the mean field phase diagrams both quantitatively 
and qualitatively. 

Our results for the 2D t-t'-V model were obtained with MATLAB using the system size 
L = 100 (i.e. 100^ lattice sites). We checked that this is large enough so that finite size 
effects are essentially negligible. However, we note that some phase boundaries are slightly 
affected by finite size effects even at this system size, as discussed in more detail below. 



3.2 2D Luttinger model 

We use HF theory as explained for the 2D t-t'-V model in the previous section. The 
reference Hamiltonian can now be written as (using a convenient matrix notation) 



Hn, = y (ct.(k),ct (k)) + ?o + gi -/i. A \ /c4k)\ 



(37) 



keA^ 

with variational parameters qo, qi and A 



The grand canonical potential corresponding to this reference Hamiltonian, evaluated 
for the Hamiltonian in f|T3|) . is denoted by Qa] see Appendix IC.2.31 for explicit formulas. It 
is important to note that this is only the antinodal contribution and that the total grand 
canonical potential of the 2D Luttinger model is 

n = na + Sa (38) 

with the energy constant Sa in fl22|) - fl26|) taking into account the contributions from the 
other fermion flavors (r, s), r = ± and s = ±, 2. Note that Qa is a function of fia (rather 
than fi), and that the filling constraint following from our general discussion of HF theory 
in Appendix [B] is 

On the other hand, the filling constraint in f l3T|) should still hold true but with ^hf replaced 
by Vt in fl551) . This implies the following consistency condition 

which must be fulfilled for the parameters /ia and £a given in Section 12.21 This identity 
is non-trivial and provides an important check of our computations. We therefore include 
details of its verification in Appendix IC.3I 
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The calculation of the grand canonical potential of the antinodal CDW phase is done 
as follows (the normal phase is treated identically). For fixed n, we make an ansatz for 
the antinodal filling Solving (fT9l) for Qo then gives us using (fT6|) . We proceed by 
minimizing the grand canonical potential f l38p with respect to the variational parameters 
go, Qi and A. This in turn gives us a specific value for the antinodal filling to be compared 
with our initial guess. We repeat the above procedure until a self-consistent solution is 
obtained for z/^. An example of a resulting curve for Q vs. fi is given in Figure |H 




Figure 4: Hartree-Fock energies vs. fi as in Figure |3] but now for the 2D Luttinger model 
at K = 0.8 and Q = 0.457r {V = 4, t = 1, t' = 0). Shown are the energies of the antinodal 
charge-density-wave (CDW) and normal (N) phases. Again there is an interval /zi < /i < /i2 
for which the CDW energy is lower, but unlike Figure [3l this energy now deviates from a 
straight line (dashed in the figure). This is due to the nodal fermions that change the total 
filling u even though the antinodal fermions are half-filled and gapped. This kind of plot 
proves the existence of a partially gapped mean field phase in the 2D Luttinger model. 

For the numerical calculation of fla, we fix the number of momenta in the antinodal 
Fourier space regions ()T0|) to 6400. Then L = S0V2/k, which is large enough so that finite 
size effects are smaller than the symbol size in our figures. Furthermore, at this system size 
we can safely replace the Riemann sums in Appendix IA.2.31 by integrals. 



4 Phase diagrams 

In this section, we present and discuss mean field phase diagrams for the 2D t-t'-V model 
and the 2D Luttinger model. 
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4.1 2D t-t'-V model 



Figures [T](a), E] and [H] give various phase diagrams of the 2D t-t'-V modeL Shown are the 
phase boundaries of the CDW phase, the N phase, and the mixed phase (indicated by 
horizontal hues) and how they depend on filhng u, couphng V and nnn hopping t'. 

As seen in Figure [T|^a), for t' = the CDW phase is only stable at half-filling u = 0.5, 
and there is a significant mixed region away from half-filling. The size of this region grows 
with increasing V. The invariance of the phase diagram under z/ — 1 — z/ is a consequence 
of particle-hole symmetry and t' = 0. We found that all phase boundaries in this figure are 
insensitive to finite size effects. 

The dependence of the phase boundaries on t' and u for = 4 is shown in Figure [51 
Due to particle-hole symmetry this phase diagram is invariant under {t', v) — )■ (— t', 1 — z/), 
and we therefore only discuss t' < 0. For < —t' < 0.15(3), the effect of t' is small and, in 
particular, the CDW phase exists only at half-filling. However, for 0.15(3) < —t' < 0.21(3) 
it is possible to dope the CDW phase on the particle (but not on the hole) side, as discussed 
in Section 13. 1[ For even larger values of —t\ the CDW phase can be doped both on the 
particle and the hole side, and the mixed region becomes smaller with increasing —t' and 
eventually vanishes. Note that the phase boundaries between the N and the mixed phases 
do not change much with —t\ for example, z/n,i = 0.70(1) and 0.69(1) for —t' = and 
0.3, respectively. Moreover, the phase boundary between the CDW and the mixed phases 
at the hole side for —t' > 0.21(3) is, to a good approximation, a straight line. The small 
wiggles of the phase boundary between the CDW and the mixed phases on the particle side 
for —t' > 0.18(3) are due to (minor) finite size effects, but all other phase boundaries are 
quite insensitive; the same is true for the phase diagrams in Figure [HI 
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Figure 5: Mean field phase diagram of the 2D t-t'-V model at zero temperature: next- 
nearest-neighbor hopping t' vs. filling z/ at t = 1 and V = 4. Shown are the charge- density- 
wave (CDW), normal (N), and mixed (horizontal lines) phases. 

The dependence of the phase boundaries on t' can also be seen in Figures M^a.) and (b) 



16 



showing the V vs. u phase diagrams for t' = —0.2 and —0.4, respectively. Note that, for 
non-zero t', there is a critical coupling value Vc > below which no CDW phase exists 
(e.g. Vc = 0.8(1) for t' = —0.2), and that there is a phase boundary between the CDW and 
N phases. Moreover, for strong coupling, the CDW phase broadens out and increasingly 
dominates the phase diagram (only visible in (b)). The phase boundary between the CDW 
and N phases is quite sensitive to finite size effect, and it is difficult to determine from our 
numerical data if it is a first- or second order phase transition. Finally, to see how the phase 
diagram evolves with t' it is instructive to compare Figures [T](a) with Figures |6]^a) and (b). 
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Figure 6: Mean field phase diagrams of the 2D t-t'-V model at zero temperature: coupling 
V vs. filling z/ at t' = —0.2 in (a) and t' = —0.4 in (b) {t = 1). Shown are the charge- 
density- wave (CDW), normal (N), and mixed (horizontal lines) phases. 



4.2 2D Luttinger model 

Apart from Figure [U^b), all our phase diagrams for the 2D Luttinger model have been 
computed for t = 1, V = 4 and t' = or —0.2. These parameter choices are partly 
motivated by our results on the 2D t-t'-V model. To be specific, for t' = —0.2, the phase 
diagram displays several, qualitatively different, features near half-filling, while V is still 
of the same order of magnitude as t; see Figure [6]J^a). Furthermore, there is an interesting 
transition near V = 4 and t' = —0.2 at which it becomes possible to both particle- and 
hole-dope the CDW phase. Other than this, there is nothing special about these parameter 
values. 

As mentioned in the introduction, a main result of this paper is that the 2D Luttinger 
model indeed has a phase in which the antinodal fermions are gapped and half filled, as 
conjectured in [T^ 2]- Figure [U^b) shows one example with fixed values of Q and k (as 
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explained later, we can in fact eliminate the parameter Q in this figure by the requirement 
Q = Qo)- Similar to the 2D t-t'-V model, we again find a CDW phase, a N phase, and a 
mixed phase in-between. However, for k < 1, the mixed phase is typically much smaller 
than for the 2D t-t'-V model; cf. Figure (H^a). 



(a) t=1 , t'=0, V=4 and k=0.7 



(b) t=1,t'=0, V=4 andK=0.7 
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(c) t=1 , f=0, V=4 and 2QM=1 .0 
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(d) t=1 , t'=0, V=4 and 2Q/jt=0.6 



Q 

. O 



0.1 0.3 0.5 0.7 
V 



0.9 



Figure 7: Mean field phase diagrams of the 2D Luttinger model at zero temperature, t = 1, 
t' = 0, = 4, and for different values of the antinodal region size k and nodal linearization 
point Q. Shown are the charge- density- wave (CDW), normal (N), and mixed (horizontal 
lines) phase boundaries of the antinodal fermions. (a) Q vs. total filling z/ at k = 0.7 (results 
for other values of k are similar), (b) Q vs. nodal Fermi surface location Qq at the four 
phase boundaries in Figure (a), (c) k vs. z/ at Q = O.Stt. (d) k vs. z/ at Q = O.Svr. Note that 
we have to restrict to k > |1 — ^| as indicated by the dashed line. Due to particle- hole 
symmetry, the corresponding result for Q = O.Ttt can be obtained as u ^ 1 — u. 

An important question is how sensitively the results depend on Q and k. We find that 
the qualitative features of the phase diagrams are robust and the quantitative dependence 
on Q is weak. However, the quantitative dependence on k is more pronounced. To be more 
specific. Figures [7l^a)-(d) give representative examples for f = 0: (a) shows how the phase 
boundaries depend on Q if k is fixed, while (c) and (d) show how they change with k, for 
fixed Q. We note that the almost vertical phase boundaries found in Figure El^a) for Q 
near 7r/2 is not special to the current choice of k and t' . Particle-hole symmetry and t' = 
imply that the phase boundaries are invariant under {Q, v) — )■ [n — Q,l — v). This and 
Q = -n /2 explain the symmetry of Figure [Tl^c). Recall that k is restricted by (fT2|) to lie in 
the range 11 — — I < k < 1. 

Figure El^b) shows Q as a function of the nodal Fermi surface location, parameterized 
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Figure 8: Mean field phase diagrams of tlie 2D Luttinger model at zero temperature, t = 1, 
V = 4 and t' = (left) respectively t' = —0.2 (right): antinodal region size k, vs. filling 
u. Shown are the charge- density- wave (CDW), normal (N), and mixed (horizontal lines) 
phases of the antinodal fermions. The parameter Q is determined such that the nodal 
points are located on the nodal Fermi surface arcs at each phase boundary, i.e. Q = Qq. 



by Qo, at the four phase boundaries in[71^a). These figures suggest that one can fix Q by 
the requirement Q = Qo using the following iterative procedure: given Q = Qn one can 
compute Qo by solving the mean field equations and then set Qn+i = Qo- We generally 
find that the sequence {Qn}n=i,2,... converges quickly independent of the starting value for 
Q. One can thus eliminate the parameter Q and obtain phase diagrams depending only on 

Figure [S] shows two examples of such phase diagrams, the left for t' = and the right 
for t' = —0.2. Comparing the left diagram with Figure El^c), one finds that the filling 
values at the four phase boundaries agree up to an error ±0.002. Since Q varies over an 
extended interval in Figure [H while it is fixed at 7r/2 in [Tl^c), this further demonstrates 
the insensitivity of the phase boundaries to changes in Q. The same feature holds true for 
the right diagram. Moreover, Figure M shows again the qualitative changes of the phase 
diagram induced by non-zero t'. For t' = —0.2, the hole side of the phase diagram is similar 
to the one for t' = 0, but on the particle side, one no longer finds a mixed region. Instead 
there is a continuous transition between the CDW and the N phase. 

Figure [H] compares the effect of varying the temperature on the phase diagrams of 
the 2D t-t'-V model (a) and the 2D Luttinger model for k = 0.8 and Q = Qo (b). For the 
2D t-t'-V model at finite temperature, it becomes possible to dope the CDW phase away 
from half-filling. Moreover, the mixed phase decreases in size as temperature is raised from 
zero, and it completely disappears at 1/(3 = 0.50(2). For larger values of there is a 
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Figure 9: Comparison between the temperature vs. filling (u) mean field phase di- 

agrams of the 2D t-t'-V model (a) and the 2D Luttinger model (b) for t = 1, t' = 
and V = 4. In (b), k = 0.8 and Q is fixed by the condition Q = Qq at the phase bound- 
aries. Shown are the charge- density- wave (CDW), normal (N), and mixed (horizontal lines) 
phases. Note the different temperature scales in the two plots. 

continuous transition between the CDW and N phase. The qualitative features of the 2D 
Luttinger model at non-zero temperature are quite similar to the 2D t-t'-V model, except 
that the CDW phase is only partially gapped and the overall temperature scale is reduced. 
For example, the mixed phase now disappears at = 0.25(2). We note that in computing 
the filling contribution from the nodal fermions, we have assumed for simplicity that there 
are only bosonic excitations from the ground state, i.e. the total number of nodal fermions 
is independent of temperature. We leave it to future work to investigate whether or not 
this is a justified assumption. 

Finally, Figure [10] compares phase diagrams of the 2D Luttinger model obtained using 
the exact antinodal band relation in fl27|) (open circles) to the ones obtained using the 
Taylor series approximated ones in ( ITTj) (full circles). Note that, for t' = (left), there 
are only small quantitative differences. For t' = —0.2 (right), however, there are larger 
deviations, in particular for u > 0.5 and k > 0.6: in the former case we find a mixed region 
between the CDW- and N phases, but in the latter case this mixed region is absent. 



4.3 Discussion 



In this section, we discuss some general features of the mean field results for the 2D 
Luttinger- and 2D t-t'-V models that we observed in our numerical computations. These 
observations are also based on phase diagrams not included in the present paper. 
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Figure 10: Comparison between the results in Figure |8] obtained with the Taylor series 
approximated band relations f|TT]) (full circles), and the ones obtained using the full band 
relations fl27|) (open circles). Note that the results agree well for t' = 0, but there are some 
differences for t' = —0.2 and n > 0.6. 

Whenever the CDW and the N phase share a phase boundary in a diagram (i.e. there 
is no mixed region in between), the CDW order parameter A goes continuously to zero 
at the boundary. Such phase boundaries are more difficult to determine numerically than 
transitions to a mixed phase. 

For t' = we only found a CDW phase in the 2D Luttinger model with Ua = 0.5, i.e. 
the antinodal fermions are half-filled. However, as for the 2D t-t'-V model, for non-zero t' 
it is possible to have a CDW phase with Ua 7^ 0.5. In fact, as a general rule of thumb, if the 
parameters V and t' are such that the CDW phase in the 2D t-t'-V model can be doped, 
then the 2D Luttinger model has a partially gapped phase with z/^ 7^ 0.5. The converse is 
not always true. We expect that the physical properties of the 2D Luttinger model in a 
gapped phase with z/^ 7^ 0.5 is qualitatively different to one with Ua = 0.5 . 

Furthermore, as exemplified in Figure [TUl when t' = the phase diagram of the 2D 
Luttinger model is hardly changed if one replaces the Taylor series approximated band 
relation in (fTTj) by the exact one in (127|) . For t' ^ and k close to 1, the results are more 
sensitive to this replacement. 

We note that it is not obvious that the phase boundaries can be fixed unambiguously 
by the requirement that Q = Qo since the Qo-value for which this occurs is, in general, 
different for the N and the CDW phase. However, we found that using the Qo-value from 
the N and the CDW phase leads to results that are very similar: for t' = the discrepancy 
is typically smaller than the symbol size in our figures, and this is also true for t' = —0.2, 
apart from the case when k is close to one and the full band relation in fl27|) is used. In 
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Figures IHlfTOl we determined the phase diagrams using Qq from the CDW phase. When 
using the Qo-values from the N phase, the CDW phase increases shghtly in size, while the 
N phase decreases shghtly in size. 

We now discuss how the size and location of the nodal arcs evolve in the left diagram 
of Figure [8] as k and v are varied. This serves as a representative example for the general 
case. Consider first fixed k and v < 0.5 {v > 0.5 is analogous). When the system is in the 
half-filled and partially gapped CDW phase, one finds Qo = As the nodal fermions 

are hole doped, Qo decreases until eventually the CDW phase becomes unstable at some 
filling u < 0.5. For even smaller values of u, the system is first in a mixed phase, followed by 
a N phase. When the system goes from the CDW-mixed phase boundary to the N-mixed 
boundary, the Qo value will first increase towards tt/2. As an example, when k = 0.8, Qq 
is 0.38(1) on the former boundary and 0.44(1) on the latter. Furthermore, when the size 
of the nodal arcs is increased (i.e. decreasing the value of k), the value of Qo on the phase 
boundaries approaches asymptotically 7r/2. 

Finally, in mean field theory, the antinodal fermions in the 2D Luttinger model behave 
very much as the fermions in the original 2D t-t'-V model scaled by a factor k^. To give 
an example, when the antinodal fermions are half-filled (z/^ = 0.5), the size of the CDW 
gap is, to a good approximation, proportional to k^. Likewise, the qualitative features of 
the temperature vs. filling phase diagrams for the 2D t-t'-V - and 2D Luttinger model are 
almost identical if the temperature scale of the latter is reduced by a factor ^ k^. 



5 Final remarks 

1. Mean field theory is a variational method and therefore not necessarily restricted to 
weakly coupled systems. For example, there exist interacting fermion Hamiltonians of 
the type ([T]) for which mean field theory is exact (examples include Hartree- and BCS- 
like models; see e.g. [TS] for details). Moreover, for many models describing electrons in 
conventional 3D metals, it is known that mean-field type approximations can give accurate 
results despite of the presence of strong Coulomb interactions. Nevertheless, for lattice 
fermion systems of Hubbard-type, standard mean field theory fails in a large part of the 
parameter regime [7] . In this paper, we demonstrated that it is possible to circumvent this 
problem by treating parts of the fermion degrees of freedom exactly using bosonization, as 
proposed in [T[|2]. It is straightforward to extend this approach to the 2D Hubbard model 

2. The results in [2] and the present paper suggest that the 2D t-t'-V model has a qual- 
itatively different behavior in different filling regimes that can be described by different 
effective Hamiltonians. For example, at half filling there is, on the mean field level, a 
fully gapped CDW phase that is adequately described by the antinodal Hamiltonian in 
ffT^ with K = 1. Upon doping the system, both the nodal- and the antinodal degrees 

^As noted in Section U7I\ and further discussed in [5], this parameter regime must be interpreted with 
some care since one has additional back-scattering interaction terms not included in the 2D Luttinger model 
when Qo — 7r/2. 
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of freedom become relevant, and the low-energy physics is governed by the 2D Luttinger 
model with k < 1. For the special case when the antinodal fermions are gapped, it is 
possible to describe the system by a pure nodal fermion model that can be solved exactly 
by bosonization [2]. For large filling {ly close to one), the low-energy physics is expected to 
be dominated by the out-fermions with r = + and s = 2. In this regime, the appropriate 
effective model describes non-interacting fermions with a band relation £:+^2(k) oc [kf + 
[2]. Corresponding statements hold for the in-fermions at small filling (z/ close to zero). 

3. A main result of this paper is that the 2D Luttinger model indeed has a partially gapped 
phase with gapless nodal fermions and gapped antinodal fermions. The possibility to obtain 
such a phase is insensitive to changes of parameter values and model details like the band 
relation. 

4. For K = 1 all degrees of freedom in the 2D Luttinger model are treated in mean field 
theory (there are no nodal fermions then). One should therefore expect that the phase 
diagram of the 2D Luttinger model for n = 1 should be qualitatively similar to the one 
of the 2D t-t'-V model. Our results show that this is indeed the case, but there are some 
quantitative differences; compare Figures and [HI Most of these differences are explained by 
the Taylor series approximation of the antinodal bands (fTTj) : cf. Figure [TOl The remaining 
discrepancy is due to the approximation of the interaction vertex mentioned at the end 
of Section 12.21 As discussed there, it is possible to improve on this approximation [2] and 
derive a refined 2D Luttinger model that gives back the 2D t-t'-V model when setting k, = 1. 

5. There exist parameter values such that the renormalized antinodal coupling constant Qa 
in (fT7|) is negative: we found, for example, ga < for t = 1, t' = —0.2, V = 10, k = 0.8, 
and 2Q/tt = 0.24. However, this and all other such cases we found are barely within the 
domain of validity of our method; recall (I12p . We therefore restricted our discussion in this 
paper to ga > 0. 

6. It is worth noting that our results do not violate the Luttinger theorem ^U\: the 
proof of this theorem assumes a standard connected Fermi surface, and it therefore requires 
modifications if there is a (partial) gap. 

7. In this paper, we only presented mean field results for the simplest consistent charge- 
density- wave ansatz (i.e. only qq, qi and mg non-zero, see Appendix [C]). Another inter- 
esting possibility is the so-called d-wave charge-density-wave (DDW) phase that has been 
suggested in the context of high-temperature superconductors [22]. It corresponds to an 
ansatz with qo, qi and 1712 non-zero (this is a special case of the general ansatz in fl73|) ). One 
can find mean field solutions for which m2 7^ gives lower energy than m2 = 0. However, 
we never obtained a DDW solution that has lower energy than the CDW solution discussed 
in the main text. The same result hold for the other mean field order parameters given in 
Appendix O 

Note added: Recently, numerical simulations of the 2D t-t'-V model using fermionic PEPS 
f2T\ produced a phase diagram that is remarkably similar to our Figure (H^a); cf. Figure 22 
in [a]. 
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A Model details 

The results in the present paper are a shght generahzation of the ones derived in |2]. In 
this appendix we explain how these generalizations are obtained. We also point out the 
(minor) differences in notation as compared to [2j. In the following, we write "(1.33)" short 
for "Equation (33) in [2]" etc. 

In [2], the lattice constant a was an important parameter when taking the partial con- 
tinuum limit in the nodal region. However, in the present paper we can without loss of 
generality set it to 1. The relation between the fermion operators used here and the ones 
in [2] is as follows, 

c(k) = f^(k), c±(k) = f V'±,o(k). (41) 
A.l 2D t-t'-V model 

The Hamiltonian in ([3]) is obtained from the one in (1.19) and (I.25)-(I.28) by inserting 
u(p) = l^M(p)/(87r^) and fHTj) . and using 

^ ^ki,k2,k3,k4C^(ki)c(k2)c'^(k3)c(k4) = ^ t;kj,k2,k3,k4c'^(ki)c^(k3)c(k4)c(k2), 

kjGBZ k^eBZ 

i.e. normal-ordering of the interaction with respect to the trivial vacuum does not generate 
any additional terms. 

We recall that the 2D t-t'-V model is invariant under the following transformation of 
parameters: 

(t, t\ V, /i, u) ^ (t, -t', 2r - /i, 1 - z/). (42) 
To see this, we make the particle-hole transformation 

c(k) ^ c"^(-k + Q), c'^(k) c(-k+ Q), t' ^ -t' , ^ ^ 2V - /j 
in dS]) and ^ and find 

N^L^-N, H ^ H + {V - iJ,)L\ (43) 
A. 2 Effective antinodal Hamiltonian 

The 2D t-t'-V Hamiltonian and the effective antinodal Hamiltonian (1131) are related 
through a series of steps which we summarize schematically as 

Hu'V ^H^ + Ha + Hna + ^0 ^ ^eff + Sn + Sq = Ha + £a- (44) 



24 



The approximate equality in (jH]) is the derivation of the 2D Luttinger model from the 
2D t-t'-V model; see Section 5 in [2]. An important feature in this derivation is normal- 
ordering with respect to a non-trivial reference state |vac). This leads to both constant shifts 
— )■ ^r,s of the chemical potentials for the different fermion flavors r = ±, s = 0, ±, 2, 
and to an additive energy constant £q. The arrow in fl44p involves integrating out the 
bosonized nodal fermions in the partition function for the 2D Luttinger model. This gives 
an additional energy constant £n and, after a local-time approximation [2], an effective 
Hamiltonian H^^ for (normal-ordered) antinodal fermions. The normal-ordering is undone 
on the right-hand side of the equality in (jB]), leading to the antinodal Hamiltonian ( |T3l) 
and the energy constant ( l22l) . 

The above-mentioned constants play an important role in the present paper, and we 
therefore give the details of how to derive them here. Note that £:r,o(k) in (1.47) and (1.6) 
is identical with £:r(k) in flTT]) (recall that k± = {ki ± k^j^l^. Moreover, Aq in (1.30) and 
( !T0|) are the same. 

A. 2.1 2D Luttinger model 

In Section 5.3 in [2j, the chemical potential /i was fixed by the condition that the one-particle 
states in the nodal regions s = ± were completely filled up to the points Qr.,± = (rQ, ±rQ). 
In the present paper the chemical potential is treated as a free parameter, and the nodal 
fermions are thus filled up to some points (rQo; i'^Qo); with determined by /x; see fH9|) 
below. With this generalization. Equations (1.59), (1.60), and (1.9) are modified by the 
replacement Q — > Qo- 

Inserting (1.53), (1.54), and (the modified) (1.59) in (1.63), we obtain for the effective 
chemical potentials for the different fermion flavors 

/^±,o = - 2Vz/ - 4t' 

;x±^2 = - 2Vz/ + 4t' T [4t + 2FC] (45) 
/^r,± = /i - 2Vz/ + 4t cos(g) + 4t' cos^(g) + WC cos(Q) 

with C in (l21]) and v in ( TT5l) . 

As argued in Section 5.3 in |2], the 2D t-i!-V Hamiltonian can be approximated by 
the normal-ordered 2D Luttinger Hamiltonian without changing the low-energy physics: 
Htvv ~ Hn + -f^a + -f^na + ^o! sec (1.62). The energy constant appearing here is 

^0 = 5^ 5^ £r,s^£int (46) 
r=± s=0,±,2 

with 

^r,s = Yl (^)MQr,s) + £.,.(k)](vac|^t,(k)^,Jk)|vac) (47) 

keA*^ 

and Sint in (1.65). The sets A*^ are the regions of the Brillouin zone corresponding to the 
(r, s)-fermions, and they are defined in (1.42). Note that e(Qr,s) + £^r,s(k) ~ e(Qr,s + k) are 
the effective band relations of the (r, s)-fermions; see (1.46), (1.47) and (1.6). The constants 
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£r,s thus correspond to the expectation values of the kinetic energy of the (r, s)-fermions 
in the reference state |vac). The constant £int arises when the interaction and the term 
proportional to — /x are normal-ordered. 

A. 2. 2 Effective antinodal model 

As shown in [2] Section 6.3, integrating out the nodal fermions amounts to replacing the 2D 
Luttinger model by the normal-ordered effective antinodal Hamiltonian: Hn + Ha + H^a — > 
i^eff+^^n- The additional constant is (in the zero temperature limit) the ground state energy 
of the nodal Hamiltonian, £n = (Hn)- It can be computed explicitly |2], but in the present 
paper we only need its /i-derivative given in ( 126|) . 

In the present paper we work with Ha in ( !T3|) without normal-ordering. To find the 
precise relation between the Hamiltonians Hcs and Ha, we use (1.37), (1.8) and (1.47) to 
write (1.91) as 

^eff = E(x)'E[^^.o(k) -H :^,,o(k)V^^,o(k): + 

' ^ . . (48) 

+2{2V - g,s) 5^(i)V+,o(p) J-,o(-p). 
p 

Inserting J±,o(p) = P±,o(p)^- Sp^oL'^u±fi, (1.46), (1.50), (SI]) and (1.59), we obtain Hes = 
Ha - S-^o - £+,0 + £u with Ha in ([13]) and the parameters Qa in ([T7D, /la = /io + Qa^af^'^, £±,q 
in dS]), and £:„ in (EH). Inserting (05]) gives (fT6D. 

To compute Qq we consider the nodal fermion branch r = s = + with the band relation 
e+,+(k) given in (1.47) and (1.6). The condition that these fermions are filled up to (Qo? Qo) 
is equivalent to £:+,+ (k) — = for k + {Q, Q) = (Qo, Qo), i-e. 

V2vFiQo -Q)- = (49) 

with vp in (1.6). From this we can compute fi in ( 1T9|) . and inserting ( IT9l) in ( l45l) we obtain 
(yUo = yU±,o) 

/xo = y2vF(go - Q) - 4t cos(Q) - 4t'[l + cos^{Q)] - 2VC cos(Q) (50) 
which generalizes (1.7). 

A. 2. 3 Antinodal energy constant £a 

As discussed above, the two Hamiltonians studied in this paper are related as Hu'v ~^ 
Ha + £a with the total energy constant £a = £n + £o — £+,2 — £-,2 + £u- Inserting ( 146|) and 
using ( H5|) we can write this as in ( l22l) with 

£Mn = £+,2 + £-,2 + 4^+,+ (51) 
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where we use £+^+ + _ + + _ = due to symmetry. Note that £±fl in (H6 

drops out due to undoing the normal-ordering of the antinodal fermions. 

In the hmit of large L, we can compute E^^sl^^ as Riemann integrals (see fl47l) ). 



{2n) 



[e(Qr.,.)+£r,s(k)], (52) 



with Sr^s ^ A*^ the set of all occupied states k in |vac). The (+, 2)-fermions are completely 
empty, and thus £+^2 = 0. The (— , 2)-fermions are completely filled, i.e. = ^-2 ■ Using 
(1.42), (1.47) and (1.6) we get 



L2 



/ ^h4(t + t') + (t + 20(fc^ + A:!)]. (53) 



The (+, +)-fermions are filled up to k + (Q, Q) = {Qo, Qo), i-e. 5*+^+ contains all k G A* 
with k+ < V2{Qo -Q). We thus obtain from (1.42) and (1.47) 



L2 



/ -—[-4tcos{Q)-4:t'cos\Q) + VFk+] (54) 

J -{k7T + 2Q-Tt)/V2 2^ J-{l-ti)TT/V2 2^ 



with vp in (1.6). Computing these integrals we obtain the result in ( l23l) . Finally, for Sint 
we insert (1.53), (1.54) and (1.59) in (1.65) and obtain by straightforward but somewhat 
lengthy computations ( l25il . 



B Hartree-Fock theory: Generalities 

For the convenience of the reader, we collect here the main facts about Hartree-Fock theory 
at non-zero temperature. A more detailed account can be found in [H [23]. 



B.l Thermal equilibrium states 

Many interacting fermion models, including the ones discussed in the present paper, are 
given by a Hamiltonian as defined in ([1]). The restriction to a finite number Af of one- 
particle quantum numbers amounts to having both short- and long distance cutoffs in the 
model. 

The thermal equilibrium state of ([1]) at inverse temperature (3 is obtained as the mini- 
mum over all density matrices p of the grand canonical potential 

nip)'^TT{pH)-^Sip) (55) 

under the particle number constraint 

(N) = -^ = No (56) 
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with "Tr" the trace in the fermion Fock space, Nq some given non-negative integer and 

^(p)^='-Tr(pln(p)) (57) 

the entropy of the state. The (unique) global minimum of Q{p) is obtained for the Gibbs 
state 



B.2 Unrestricted Hartree-Fock theory 

HF theory at non-zero temperature amounts to restricting the search for the minimum of 
fi(p) to the set of all HF Gibbs states defined as followsjll 

g— /3i^HF 

= Tr(e-/^^HP) (^9) 

with 

= Yl Kdci (60) 

kl 

a Hamiltonian for non-interacting fermions and hki the matrix elements of a self-adjoint 
TV X A/" matrix h. We refer to h as variational one-particle Hamiltonian. 

Note that we do not consider the most general definition of a HF state; one could also 
allow for states for which U{1) gauge invariance is broken (see the remark in lB.Sp . or states 
that are combinations of a pure state and a HF Gibbs state (see |1] for details). 

One can compute the grand canonical potential ^hf = ^(phf) for such HF Gibbs states 
in terms of the eigenvalues ex and corresponding orthonormal eigenvectors fx of h as follows, 



klmn 

1 (3ex 

' P 2^\QPe^ + 1 



kl klmn / r^-, \ 

(61) 



B ^\ e^<=> -H 1 ^ 



with 



dc 

A 



=T^^(/a).(/a). (62) 



and where {fx)j are the components of the eigenvector fx- Moreover, the particle number 
constraint is 

In unrestricted HF theory one does not make any restriction on the variational one- 
particle Hamiltonian h, and thus one has, in principle, N"^ real variational parameters. 
However, the number of variational parameters can often be reduced by the following: 



^AU summation indices in this section go over 1,2,... ,N . 
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Proposition: Local extrema of ft{p) in the set of HF Gibbs states are acquired for a 
Hamiltonian 

kl klmn 

(64) 

with (c^C;) = i-&- when (■) is the expectation value in the state pup corresponding to 
»• 

The proof is outlined at ttie end of tliis section. Tlie corresponding result for zero 
temperature can be found in 

It follows that one can restrict the search of the minimum to variational one-particle 
Hamiltonians of the form Hm = t^i — fiSki + wm, with 

'^kl ~ i'^knlm + "^nkml ~ '^knml ~ '^nklrnjlmn (65) 

mn 

for some possible (we will come back to this in Section [6.31 below). For the models of 
interest to us, this fact allows to reduce the variational parameters in HF theory consid- 
erably. For example, for the 2D t-t'-V model on a square lattice with sites one would 
naively expect variational parameters in unrestricted HF theory, but using ( 165|) this 
number can be reduced to 5L^. 

Note that together with (c|,q) = jik gives a self-consistent system of equations 
that is often used in practical implementations of HF theory. However, when restricting 
HF theory it is important to check that a self-consistent solution found in this way is indeed 
an absolute minimum of f l6ip in the set of considered states. Moreover, in restricted HF 
theory it is important to work with the grand canonical ensemble (i.e. fix /i and not the 
fermion number) even at zero temperature, as explained in the main text. 

The 7jfc are the components of a self-adjoint matrix 7 = (e'^^' + l)^^ called a one-particle 
density matrix, and this matrix completely specifies the corresponding Hartree-Fock Gibbs 
state. Moreover, the results in fj63|) and the second line in fl6T|) are equivalent to tr(7) = A^'o 
and 

^S{puf) = tr(/i7 + (1//3) ln(l + e"'^'^)) = -ltr(7ln(7) + (1 - 7) ln(l - 7)) (66) 

with "tr" the J\f x JV matrix trace. Finally, in the zero-temperature limit /3 — )■ 00, the 
one-particle density matrix fl62|) becomes ■yjk = J2x^i~^)^)if^)jif>dk] see also fl30|) . 

Proof of the Proposition: To find the variation of the grand canonical potential we use (1611) 
and (IMj) . and we regard 7 rather than h as variational parameter: Qhf = ^hf(7)- We 
compute the variation 5f2 = ^f2HF(7 + 557)1^=0 and obtain 

SriuF = ^(4« - pSki + Wki - hki)5'yik (67) 

kl 

with hki the matrix elements oi h = — (1//3) ln[7(l — 7)"^] and Wki in ( l65|) : the first three 
terms are obvious, and the last is obtained from the variation of the entropy term in f l66l) . 
This implies the result. □ 
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B.3 Restricted Hartree-Fock theory 



The Hamiltonian ([T]) often has a large symmetry group, and assuming that the HF Gibbs 
state is invariant under (a subgroup of) these symmetries allows to reduce the number of 
variational parameters even further. For example, the 2D t-t'-V model is invariant under 
translations, parity, and (discrete) rotations. If one restricts to HF Gibbs states that are 
invariant under all these transformations, the number of variational parameters can be 
reduced to just two; see Section [UJI However, as discussed in the main text, symmetries of 
the Hamiltonian are often broken, and it is therefore important to also allow for HF states 
where (some of) the symmetries are broken. 

In practice, one assumes that the one-particle density matrix commutes with some 
suitably chosen subgroup of the original symmetry group. This leads to restrictions on 
the matrix elements jki, and thus, through fl65l) . on the matrix elements of the variational 
one-particle Hamiltonian h. The latter is true since the proposition stated in Appendix IB. 21 
can be generalized to restricted HF-theory (we plan to present details elsewhere). 

Remark: We note that the Hamiltonian ([1]) is invariant under the gauge transformation 

c, ^ e-c, (68) 

for real parameters a. There exists a natural extension of HF theory allowing also for 
states for which this gauge symmetry is broken. This extension is relevant for models with 
attractive interactions, and its physical interpretation is e.g. superconductivity!^ However, 
since the interaction of the 2D t-t'-V model is purely repulsive, superconducting HF states 
cannot occur; see |1] for proof. For the effective antinodal model in ( IT3l) there do exist 
parameter regimes where the coupling constant Qa is negative and thus a superconducting 
HF state is possible. However, as discussed in Section [5l Remark 5, we found that this 
parameter regime is very small. We therefore ignore superconducting HF states throughout 
this paper. 

C Mean field theory 

We give here some additional details on our mean field treatments of the 2D t-t'-V - and 
2D Luttinger models. 

C.l 2D t-t'-V model 
C.1.1 Symmetries 

The 2D t-t'-V Hamiltonian in is invariant under the symmetry group generated by the 
following transformations 

Tj : c(k) ^ e^^^c(k), V : c(k) -> c(-k), 7^ : c(k) = c{ku fcs) c{k2, -h) (69) 

^There are other physical interpretations in the context of nuclear- and elementary particle physics. 
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for j = 1,2. The transformations 7i and 72, V and TZ correspond to translations by the 
lattice vectors ei = (1,0) and 62 = (0,1), parity transformation, and rotation by 7r/2, 
respectively. 



C.1.2 Restricted Hartree-Fock theory 



We consider HF theory restricted to variational states for which the symmetry is broken 
down to translations by two sites, i.e. the restricted symmetry group is generated by (7i)^ 
and (72)^- In such a state, one has for the one-particle density matrix 

7(k', k) = (ct(k)c(k')) = e(k) 5k,k' + e(k) Vk'+Q (70) 

with Q = (vr, vr) and some functions 6 and 6 satisfying 



0(k) = e(k) = 0(k + 2Q), 0(k) = e(k+Q) = e(k + 2Q). (71) 

Using the fact about HF theory stated in fl64l) ff we can restrict ourselves to HF potentials 
of the following form, 

w(ki, k2) = -jj ('i^ki,k2,k3,k4 + 'yk3,k4,ki,k2 - ^^ki,k4,k3,k2 " "Wka .ka ,ki ,k4 ) 7 (^4 , ks). 

k3,k4eBZ 

Inserting flTOj) and ([S])-©, and using u{k — k') = ^ J2'j=i with 

Wi,2 = cos(A;i) ± cos(A;2), ^3,4 = sin(/ci) ± sin(A;2) (72) 

we obtain 

w(k, k') = ^go + ^ QjUjiX)^ 4,k' + ^"^0 + (^k.k'+Q (73) 



with 



and 



V V ^ 

go = 2VnQ, mo = -2l^no, qj = -^^j, rrij = -^nj (74) 



^o = ^E[Q(k) + e(k + Q)], no = -^$^Ree(k) 

k k 

= ^ E[®(k) - e(k + Q)]«,(k), n, = ^Y. ime(k)«,(k) 

k k 

for j = 1, 2, 3, 4, with all sums over BZ; we used fITT]) and Uj(k + Q) = —Uj(k). 



(75) 



31 



C.1.3 Mean field equations 



To find G and G we use Fourier transformation and write tlie reference Hamiltonian in fl2^ 
and ( 132|) in matrix form as follows, 

iJHF= E (ct(k),ct(k+Q))Mk)(^^"^^)^^ 

keBZi/2 ^ 

with 

Mk) = 

and 



a+(k ) 6(k) 
6(k) a_(k) 



a+(k) = e(k)+go-Ai+E'?jMj(k), a_(k) = a+(k+Q), 6(k) = mo+i ^ mjMj(k), (76) 

J=l 3=1 

where the sum is restricted to half of the Brillouin zone, BZi/2 = {k e BZ : ki> 0}. We 
now can compute 7(k) = f{h{k)) for /(x) = l/( + 1) using the following general result, 

..... ^^ _ /(e+(k)) + /(e_(k)) ^ /(e+(k))-/(e_(k)) / ai(k) 6(k) 
2 ^+ WiX) \W) -«i(k) 

with 

e±(k) = ao(k) ± iy(k), iy(k) = v/ai(k)2+|6(k)P, ao,i(k) = i[a+(k) ± a_(k)]. (77) 

The e±(k) are the eigenvalues of the matrix h{\i) and equal to the effective mean field band 
relations. From this and 

' G(k) Q(k) 



7(k) 
we obtain 



e(k) e(k + Q) 



BW4(.^. ' ^'-^^ ' 



2 \^ ^ I g/3e_(k) j^ij 2W{\i) V e^'^+('') + 1 e^^-W + 1 

e(i:) ^(^) f ^ ^ 

^ ^ 2iy(k) V e'^^+C') + 1 e/5--M + 1 
By straightforward computations we obtain from (!6T!) 

with 



(78) 



4 / 1 ^ 

j=o V j=i 



(79) 



-- = -i:^E^l-(l+--'^^^'^0' (80) 

r=± k ^ 
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rij and hj in (175|) . and e±(k) in (1771) . Mean field theory amounts to minimizing Ohf in 
(!79|) with respect to the 10 variational parameters qj and m^, j = 0,1,2,3,4. Note that, 
by construction, the saddle point equations dVtnv/dqj = dfluF/dnij = are identical with 
the mean field equations in (1711) - (I75I) . 

We were only able to find absolute minima of the HF grand canonical potential with 

q2 = qs = q4 = m2 = ms = 1714 = 0. (81) 

It is interesting to note that flHTl) is equivalent to parity- and rotation invariance of the HF 
state: for HF states invariant under V and TZ one has 

e(k) = e(-k) = 0(A;2, -A;i) 

and similarly for 0. This implies that rij and hj in fl75|) vanish for j = 2,3,4, and thus 
qj = rrij for j = 2,3,4. Moreover, if this is true, then the HF equations qi = —Vni/2 and 
nil = —Vhi/2 are in contradiction unless mi = 0. Note that rrtj = for j = 1,2,3,4 is 
equivalent to 

e(k) = 0(k). 

We therefore restrict our discussion in the main text to the simplified ansatz in ( 132|) with 
only three variational parameters qo, qi and mo = A. 

C.1.4 Numerical details 

We solved the mean field equations given above using MATLAB. We found that it is nu- 
merically difficult to minimize the grand canonical potential ^hf directly, but, similarly 
as for the Hubbard model [1], one can construct an auxiliary extremization problem for 
a related potential. This latter auxiliary potential has extremal^ that coincide with those 
of f^HF (in general, the minima of f^HF correspond to saddle-points of the auxiliary poten- 
tial). Furthermore, the values of these two potentials are equal at corresponding extrema. 
This auxiliary potential is numerically well-behaved and it therefore allows to reduce the 
computational time considerably (we plan to present details on this elsewhere). 

C.2 2D Luttinger model 
C.2.1 Symmetries 

While the effective Hamiltonian in f lT3|) naturally inherits the following parity- and rotation 
symmetries from the 2D t-t'-V model, 

P: c±(k) ^c±(-k), TZ: c±{k) =c±{h,k2) -^c^{k2,-h), (82) 

its translation invariance is 

% : c±(k)^ e'^-i'c±(k) (83) 
^By "extrema" we mean here points for which all partial derivatives are zero. 
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for all translation vectors x = (xi,X2) with components Xi,X2 G M. This enhancement of 
translational symmetry is due to the continuum limit taken when deriving the 2D Luttinger 
model, and it makes it impossible to distinguish between translations by an even and odd 
number of lattice sites. We thus have the following symmetry transformations of H in f ll3p . 

C„ : c±(k) ^ e±'"c±(k) (84) 

with real a. We use particle physics terminology and refer to as chiral transformation. 

C.2.2 Restricted Hartree-Fock theory 

It is known that chiral symmetry breaking in a continuum model corresponds to breaking 
the translation symmetry from one to two sites in the corresponding lattice model; see e.g. 
[21]. We thus restrict HF theory to all states invariant under all translations 7^ but not 
under chiral gauge transformations C^, i.e. we make the following ansatz for the one-particle 
density matrix 

7.v(k', k) = (c|;(k)c,,(k')) = <5k,k'e,,.(k) (85) 
for r, r' = ± and some functions Q^r' obeying 

e^w(k) = e^/^(k). (86) 

We insert (185|) in fl65|) and obtain 

H^^ = Ho + ^Y.ll ct(k)c_,,(k)rr'e..'(k') 

k,k' r,r'=it 

which is equivalent to (!37|) with 

go = 9ano, qi = -gani, A = -251^^0 (87) 

and 

no,i = ^J2^e^+i^)±e-m, no 

k 

with all k-sums over A*. 

C.2.3 Mean field equations 

The reference Hamiltonian in f l37|) has exactly the same form as the one treated in Sec- 
tion IC.1.3t and we therefore obtain the same results as in f l77|) and f lTHj) but with 

= A. (89) 

nj - A\ho\^) + £a (90) 
with chf as in (15U|) . no,i and hi in ( ISSj) and e±(k) in (I77|l and (15^ . 




a±(k) = e±(k) + qo ± qi - fia, b(k) 
Moreover, as in Section IC.1.31 

-jY = euF - qono - qiUi - Auq - Auq + ^ V^o ~ 
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C.3 Consistency check 



We now show that the condition in fHU]) holds true. Using fIFI]) . (jSSl) and flMl) we get 
^^^1^ = [-4t cos(Q) - 4t'cos^(g) + vpV2{Qo - Q)]^{1 - «:)^ 

= [/i - 2W + 2l^Ccos(Q)]|-(z/ - K^l^a) 
where we used ( HUil . and (IT^ . Equation imphes 



a/i \ a/i Ojj, J Oji oil oji 

where we used ( 1TB]1 . We use and compute, recalhng (!?!]) . 

^ = _^ _ _ 2W)|^^ - 2\/C^ = -i. - (/X - 2Vu)^ - 2VCcos{Q)^{iy - n'u,). 
oji oil on oil ail 

We finally use ( 126|) . Recalling ( 122|) and adding the results we find that many terms cancel 
and that the remaining terms add up to the r.h.s. of (l40ll . 
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